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Equalizers 
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Two basic problems in the equalization of data channels using pulse- 
amplitude modulation are considered. The first of these is to determine just 
what pulses can be equalized and what the general equalization solution is 
when using a transversal filter. The second problem is to determine if a 
simple iterative search routine will converge to a solution if its exists. 

The unequalized channel impulse response is represented by a polynomial 
whose coefficients are the sample values of the impulse response. If no 
roots of this polynomial lie on the unit circle, the channel can be equalized. 
The transversal filter which equalizes the pulse has tapiveight values given 
by weighted sums of powers of the polynomial roots. 

Various necessary and sufficient conditions for iterative convergence 
are developed. Iterative convergence can be guaranteed if the proper linear 
weighting of the output sample errors is used in adjusting the tap-weights. 

I. INTRODUCTION 

This paper is concerned with certain aspects of the automatic equal- 
ization of low-noise, linear channels which are to be used for multi- 
level pulse-amplitude modulated (PAM) signals. The purpose of the 
equalizer is to compensate for the channel transfer characteristics in 
such a way that the over-all impulse response of the channel is a 
Nyquist-I type of pulse, 1 - ■• 3 that is, as is illustrated in Fig. 1, a pulse 
with a central peak and uniformly spaced zeros with period T. If 
such an impulse response is achieved, a sequence of amplitude modu- 
lated impulses with period T can be transmitted and the sequence of 
amplitudes can be recovered at the receiver by simply sampling in 
synchronism with period T. 

Certain obvious questions such as how to achieve synchronism and 



* University of Washington, Seattle. 

1775 



1776 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1968 

^-"-CENTRAL SAMPLE VALUE 




SAMPLE TIMES WITH SPACING T 



Fig. 1 — A Nyquist-I pulse with central sample and zeros at periodic sample 
times. 

the effects of sampling jitter, nonlinear distortion, and additive noise, 
although of great importance, are neglected in this presentation in 
order to concentrate on the methods of adjusting the equalizer. Thus, 
we assume a perfectly synchronized, noiseless, linear channel with 
an ideal sampler. Some further constraints which should simplify and 
clarify the presentation are as follows. We specify that the equalizer 
is to be a transversal (tapped delay-line) filter with tap weights {« fc } 
which can be adjusted. If the input to this filter, illustrated in Fig. 2, 
is /? (t) , then the output, y (t) , is 



y(t) = Z a t p(t + jT). 



(1) 



Notice the tap-weight numbering convention and the treatment of 
the delay-line as being composed of negative as well as positive delay. 
These conventions will simplify the notation in future derivations. 
Although the equalizer may be placed at many points within the 
communication system, for convenience we will consider it to be the 
final component other than the final sampler. Thus, the objective is 



/S(t) 




•ret) 



Fig. 2 — A transversal filter with input pit), tap weights a } , and output yit). 
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to adjust the tap weights so that the output, y{t), is a Nyquist-I 
pulse when the input, p(t), is the impulse response of the system be- 
fore equalization. 

As a final simplification we assume that the tap-weight adjust- 
ment is to be carried out when the impulse response is available at 
the input. This "training period" assumption avoids the added com- 
plexities of extracting channel characteristics when data are being 
transmitted. 4 However, the properties to be developed can be readily 
extended to continuously adaptive equalization. 

II. TIMING AND THE ROOTS OF THE IMPULSE RESPONSE 

Since the output of the equalizer is to be sampled, we need con- 
sider only its effect on the output sample values {y k }. Ideally, with 
an impulse applied to the channel input, one of the output samples, 
say y , should have unit amplitude while all the others are zero. The 
tap weights are to be adjusted in order to approach this goal. An ad- 
ditional parameter which will affect the equalization is the timing. 

Let the sample values of the impulse response at the equalizer 
input be the set {/?,} where 

fit = MT) (2) 

and 

0(0 = h(t + r). (3) 

Equation (3) is to indicate that the sampling times are arbitrary. 
This is, if h(t) is the channel response to an impulse applied at t = 
0, then the sample set {/?,} is a function of the factor t. Notice that 
our assumption of perfect synchronization means that the periodicity 
factor, T, in equation (2) is the proper value. But it does not imply 
that t is prescribed. We shall see that the operation of the equalizer 
depends very strongly upon the value of t. 

This impulse response is to be equalized by the transversal filter 
with tap weights {«,,}. The transversal filter output sample set is 
{y ; J where 

7* = £«,-&+< • (4) 

The sample set {/?,} will be considered finite in extent, that is, 

0, = for j < -m and j > M. (5) 

This is a reasonable approximation for any actual channel. 
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This sample set may be treated as the coefficients of a polynomial 
B(z). 

B(z) = p u z"* m + Pm-iz"*'- 1 + ■■■ 

+ fa m + ■•• +^_ m + 1 2 + ^_ m . (6) 

In factored form B (z) may be written as 

B(z) = P„(z - 6i)(z - d 2 ) • • • (z - 6:){z - 0,)(* - 2 ) • • • (z - 0n) (7) 

where the roots inside the unit-circle are denoted by the 9 values and 
the roots outside by the «£ values. 

I e | < 1 < | <t> | (8) 

/ + fi = M + m. (9) 

As we have noted, the sample set {fi/} is a function of the factor t. 
Thus, for any particular channel, the roots of equation (7) will 
wander as t is varied. Each root will wander on some cyclic path 
which has a period T. That is 

d,(r + nT) = 6,(r) (10) 

This is illustrated in Fig. 3 where a pulse shape is shown and in Fig. 4 
where the root loci are shown for variations in t. Notice that as t in- 
creases from to T, at least one of the roots, regardless of the pulse 
shape, will cross the unit circle. 

Perhaps the periodicity of the roots can be better understood if 
the sampling is thought of as multiplication by a comb of impulses 
with spacing T. Each impulse has associated with it a power of z. For 
example, in equation (6) , we see that the impulse yielding the earliest 
nonzero sample (/?-,„) is associated with the zero power of z, the next 
impulse yielding /?_,„ + i is associated with the first power, and so on. 
As the comb is moved relative to the pulse, f3(t) , the impulses produce 
different samples and when moved a whole period T, the comb will 
reproduce the original samples again. However, each sample would 
be paired with a one-higher or one-lower power of z than previously. 
Thus, for example, if the comb were shifted by T so that the powers 
of z were one higher, the factorization of equation (7) would be ob- 
tained with the same roots except for an additional root at z = 
since the original polynomial is multiplied by the first power of z. 
As the comb is shifted along, the additional root, which must even- 
tually go to z = 0, comes in from z = oo and at some particular shift 
crosses the unit circle. 
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Fig. 3 — A pulse and four nonzero sample positions which yield a polynomial 
with roots marked 1 in Fig. 4. Nine additional sets of roots are obtained by mov- 
ing the positions above to the right in increments of T/IO. 
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REAL AXIS- 



Fig. 4 — Root loci for pulse of Fig. 3. Scale inside circle is magnified four times. 
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III. IDEAL TAP-WEICHT ADJUSTMENT 



Let us now consider the problem of forcing the y k to zero. The 
desired conditions are 



and 



7 



7* 



o = £ OfPt = 1 



= S «,-ft + * = k * 0. 



(ID 



(12) 



Equations (11) and (12) in the form of a matrix array is as follows. 
We shall assume temporarily that the number of taps is infinite. 



0/3,, ft 0, O 0-i /3- a 

0/3,,, • • • ft, ft ft ft ft., 0_ 2 

ft, • • • 03 ft ft ft 0-, 0- 2 • 

0^, ••• 03 & ft ft, (8-, 0-2 



0-„, o 

• • 0- m • • 

• • • 0- m 

(8— 
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a_2 







a-a 














(13) 



Now let us consider a typical equation of the form, 

M-k 

£ otjPt+i = 



(14) 



im-m-t 

where we will take either the case where k < —m or k > M . If these 
conditions on A; are satisfied, then equation (14) involves a values with 
only positive subscripts or only negative subscripts. Equation (14) is 
a linear homogeneous difference equation in the variable aj. Such 
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equations have solutions of the form 

(15) 
Substitution of equation (15) into equation (14) yields 

PmZ + Pu-iZ + • • • + j8i2 

+ PoZ-" + fi-.Z- 1 -" + • • • + P- m Z- m - k 

= z- m - k [0„z A,+m + fi u - lZ M+m - 1 + • • • + fa m + /3-J 

= 2 - m - A B( 2 ) = 0. (16) 

Equation (16), and hence equation (14), will be satisfied for z 
equal to one of the roots of the impulse response polynomial B(z). 
And since we have a linear homogeneous equation, any linear com- 
bination of solutions is a solution. Thus, the general solution is 

a,- m CM + C 2 d' 2 + ■ ■ ■ + Cj0i + D«t>[ + •• ■ + D^ (17) 

The constants, C±, Cz, . . . , Ci, D lt D 2 , . . . , Do, are arbitrary and 
will be adjusted to meet the boundary conditions. One boundary con- 
dition which should be imposed is the following. Eventually, we must 
approximate this infinite delay line with one of finite extent. This 
truncation should throw away only taps of small magnitude, and thus 

the tap weights should decrease in magnitude away from the center 
tap. Consequently, we demand that a } — > as h'| -> oo. Thus, 

C, = C 2 = • • • = Ci = for j negative. 

D, = D 3 = ■ ■ ■ = D a = for j positive. (18) 

In effect, we have two solutions; one for taps with negative sub- 
scripts and one for taps with positive subscripts. 

a,- = C x 6\ + C 2 6[ + • • • + CM for j > o (1Q) 

a,- = Drfl + Drfl + • • • + Drfi for j < 0. 

The region of overlap in equation (13), that is, the region where 
the equations involve tap weights with both positive and negative 
subscripts, will determine the arbitrary constants. To illustrate this, 
let us consider the impulse response of Fig. 2. For one set (No. 8) of 
samples, the samples are —2, 3, 11, —6 and the roots are 0.5, 3, —2. 
Thus 



a,- = C,(0.5)' for j > 

a . = D,(3)' + D 2 (-2)'" for j < 0. 



(20) 
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Now, C lt D lt and D 2 may be found by using equation (20) in equa- 
tion (13). However, first we must assign an origin in our pulse. That 
is, we must decide whether —2, or 3, or 11, or —6 is to be called /? - 
Solutions exist for each of these possibilities, but only one of these has 
certain desirable properties which we will discuss later. In this case, 
p should be the sample of magnitude 11. With this assignment, equa- 
tion (13) may be written 



-2[C,(.5)-'l + 3[C,(.5) 3 ] 4- 11[C,(.5) 2 1 - 6[C,(.5)] = (21a) 

-2[C,(.5n + 3[C,(.5) 2 ] + 11[C,(.5)] - 6*o = (21b) 
-2[C,(.5) 2 ] + 3[C,(.5)] + ll« n - 6[D,(3)-' + D„(-2)-'] = 1 (21c) 
-2[C,(.5)] + 3a„ + llf/),^)-' + D 2 (-2)"'] 

- 6f£>,(3)- 2 + D,{-2T 2 } = (21d) 
-2«„ + 3[/),(3r' + D 2 (-2)- 1 ] + ll[D t (Sy 2 + D,(-2y 2 } 

- 6[D>(3)- 3 + D 2 (-2)- :! ] = (21e) 
-2f/) 1 (3)- , + A,(-2)- 1 ] + 3[A(3)" 2 + D 2 (-2y 2 } 

+ 11[A(3)- 8 + D 2 (-2)- 3 } - 6[D,(3)- 4 + D a (-2)-] - (21f) 



Equations (21a) and (21f) and all others above and below these two 
are automatically satisfied for any choice of Ci, D± t D 2 . Consider 
equation (21b). In order for it to be satisfied, a must be equal to C x 
(0.5)°. Similarly, for equation (21e) to be satisfied, a© must equal 
D 1 (3)° + 7J 2 (-2)°, 



Cj = D L 4- D 2 = a,, 



(22) 



and equation (21) yields the following values for the C and D con- 
stants. 



D, = 2/95 
D 2 = 5/95 
C\ = 7/95 



& = 11. 



(23) 
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If we had chosen to let /? be the sample of magnitude 3, we would 
obtain 



"ft, = 3. (24) 



D, = — 10/95] 
D, = 15/95 

C, - 5/95. 
If the sample of magnitude 2 is made /?«, then 

D, - 9/251 

D a = 4/25L = -2. (25) 

C, = 1/50. 

In these last two examples, the solutions may be considered inferior 
because the tap weights away from the center tap will be larger than 
in the first example. However, stronger objections to the last two 
choices for /? will be raised shortly. 

A point of considerable interest is apparent in the development 
above. If any of the roots fall on the unit circle, then no solutions 
exist in which the tap weights decay in magnitude in both directions 
away from the center tap. 

IV. TRUNCATION EFFECTS 

In any practical equalizer, the number of taps available is not 
infinite. Thus, we must investigate the effects of limiting the number 
of taps to some reasonable finite value. For example, let us suppose 
that we have N + n + 1 taps. 

a, - for \ J <: ~ n - : r,i 



j> N 

Let us consider two schemes for setting the truncted tap weights. 
A more or less obvious way is simply to take the infinite solution 
[for example, equation (23)] for all available taps. This can be rep- 
resented in matrix form : 

[B*]M = [7]. (27) 

The matrix [B x ] is the infinite matrix of (3 values shown in equation 
(13) . The truncated tap set is [a T ] , 



1784 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1968 



[a T ] = 






«,v 



(28) 



and [y] is the matrix of output sample values. Notice that [a T ] may 
be written as 

M = [«.] - M (29) 

where [a*,] is the infinite set of tap weights which give us the desired 
output and 



OtN+2 





• (30) 



a-n-i 

tt-n-2 



EQUALIZER CONVERGENCE 1785 

Thus, equation (27) becomes 

[BJ\a r ] = [B.kl - WJM = M - M- (31) 

The desired output [/*], all zero samples except the central one with 
unit magnitude, is given by [B w ] [a x ] while the error set [c] is given 
by [£«][«,]. 

It is possible to calculate this error set exactly, but an approxi- 
mate bound should suffice. Assume that N and n are large so that the 
a values in [a t ] depend only on the roots which are the nearest to the 
unit circle. Let the magnitude of nearest inside and outside roots be 
6 V and <f> u . Then approximately 

a, « < fOT ' > » ■ (32) 

<j>' u for j < —n 

Consequently the largest component of any error term in [«] will be 
about 

Max comp £g aax • C u U V+ ' S 0™* Dj>Z n ~ ' (33) 

and each error term is the sum of less than N + n + 1 components. 
Thus, an upper bound on the error terms is 

I e,,,„x I < ! (A r + n)CX +1 |S|(i^+n) AA~"~'| (34) 

which vanishes as N and n are made very large. 

V. A SECOND METHOD OF TRUNCATION 

This second method offers no improvement in ultimate equaliza- 
tion over the method just discussed. However, it does lend itself to 
iterative adjustment techniques whereas the first method tacitly as- 
sumes a computation which provides the proper infinite solution to 
begin with. In this second method we require that all the output sam- 
ples (excluding y ) corresponding to the N + n + 1 taps, that is, 
jn, yjT-i, ■ ■ ■ , V2, yii y-i, y-2, • • • , y-n+i, y-»» be zero. This criterion 
may be called the Lucky criterion since it is the one R. W. Lucky 
has used in his work. 5 

What docs this criterion mean in terms of the solutions (powers of 
impulse response roots) we discussed for the infinite tap case? We are 
essentially constraining our system further by another set of boundary 
conditions. We will call these boundaries the positive boundary at ax 
and the negative boundary at a_„, in addition to the central boundary 
around a where we have already discussed satisfying boundary con- 
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ditions (those specifying d , C 2 , ■ ■ ■ , Di , D 2 ,...). As an example, 
let us consider the equations at the positive boundary. 

/W + ft,**-, + + ft m a„_ m = 7* =0 (35a) 

fta A . + fta jV -, + • • • + /3_ m aAr- (m + ,) = 7at-i = (35b) 

Pur-iOtN + (J M - 2 a N - l + ••• + /3_,„a^_ (m + w -,, = 7jv-<«-d = (35c) 

0MOCN + ftz-^V-1 + + /3-m«AT- (m + A/) = JN-Af = (35d) 

Equation (35d) can be satisfied with the solutions determined by 
the central boundary, that is 

a,- m del + C 2 d' 2 + • • • + c,e', . (36) 

However, this solution will not satisfy the M equations above equa- 
tion (35d) since the complete set of /? values is missing in these equa- 
tions. The exponentially growing solutions which were discarded 
earlier must be used now. Thus, 

a, = d0{ + • • • + CM + cM + cX + • • • + c^i for j > 

(37) 

<*,• = Drfl + • ■ • + Dati + cl { 6\ + d 2 $i -f • • ■ + d,6\ for j < 

(38) 

must be used in order to satisfy all the boundary equations. If O = M, 
then the lower case c's in equation (37) provide just enough con- 
stants to satisfy the M equations of equations (35a through c). Fur- 
thermore, J will equal m and the / lower case d's will provide just 
enough constants to satisfy the m boundary equations at the nega- 
tive boundary. 

To illustrate the preceding discussion, let us return to the specific 
example discussed previously. We use the results of equation (23) 
with 

ft - -2, ft = 3, ft - 11, ft-i = -6. (39) 

The positive boundary equations are 

HICCJO* + CiW + c 2 (-2) N ] 

- etcos)"- 1 + ^3)"- + c a (-2) J " 1 ] = o (4Q) 

3{Cr(.5) N + Cl (3)" + c 2 (-2f] + U[C 1 (.5) U - 1 + c 2 (-2) N ' 1 ] 

- 6[Ct(.S) N - a + Cl (3) N ~ 2 + c 2 (-2) N ~ 2 ] = 
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3 ,v "'(33 - 6) -2" -1 (-22 - 6) 

_3*~'(27 + 33-6) -2' v_2 (12 - 22 - 6) J 

.5 W (11 - 12)d 
_.5*(3 + 22 - 24)d. 



(41) 



If we assume that N is large enough so that the value of the C 
and D variables are unaffected by inclusion of the other roots, then 
these C and D values may be used to determine the c and d values. 
With d = 7/95, equation (41) becomes 



9 14 
.6 -4 






-e 

e 



(42) 



where 



Thus, 



3'V, 



cf 2 = (-2)% e = (r^)(-5)' v - (43) 



c\ = e/12 



c, = 



M- 5V 

(3) J 



c ' 2 = _e/8 



^ r 



(44) 



(-2)" 



Similarly, di can be found by the single boundary equation at the 
negative boundary. 



rfi = - 



Kclk" + 9 (<§>- 2r ' 



12 



(.5)' 



(45) 



The results worked out above can be roughly represented graphi- 
cally as in Fig. 5 where the magnitudes of the roots to the tap-number 
power are illustrated. This figure shows what will be called a "good" 
solution. That is, the decaying solutions predominate with the grow- 
ing solutions contributing only a small amount at the positive and 
negative boundaries. The residual errors, that is, the y k values for fc > 
N and k < —n, will be of the same order of magnitude as those of 
the first truncation method. 

In order to have a good solution as demonstrated above, the proper 
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Fig. 5 — Solution behavior for truncated equalizer. 

number of arbitrary constants must be available to satisfy the posi- 
tive and negative boundary equations. The number of equations at 
the positive boundary is M, and the number of available constants 
(outside roots) is O; similarly, there are I constants (inside roots) 
available for the m equations at the negative boundary. Thus, the 
necessary conditions for a good solution are 



no. outside roots = fi = M = no. of samples following /3 
no. inside roots = I = m = no. of samples preceding /3 



(46) 



We may use the same simple pulse to demonstrate a "bad" solu- 
tion. If we go back to equation (24), we see that an infinite (non- 
truncated) solution exists for the situation in which we decided to let 
Pa = 3. The number of inside and outside roots remain the same in 
this case, but M and m are both changed. Now M = 1 and m = 2, 
and equations (46) are no longer satisfied. All the boundary condi- 
tions can still be satisfied, but not in such a simple manner. That is 
to say, in the preceding example, satisfying equations (46) , a sepa- 
ration of solutions is possible. The central boundary specifies the 
values of the upper-case constants, then the lower-case constants are 
set to compensate for the truncation effect at the positive and nega- 
tive boundaries. Since the necessary compensation is small, and since 
the effect of the lower-case constants dies out towards the central 
boundary, only minor or negligible corrections to the upper-case con- 
stants are necessary to keep the central boundary conditions satisfied. 

Now let us consider what takes place when this step-by-step solu- 
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tion is attempted when equations (46) are not satisfied. Suppose &, 
Dj , and D 2 are set to satisfy the central boundary according to equa- 
tion (24). At the positive boundary, only one of the lower-case c's is 
needed so no difficulty arises and the positive tap solutions will, to this 
point, be little different from that illustrated in Fig. 5. However, at 
the negative boundary, the single available arbitrary constant d t is 
not enough to satisfy the two boundary equations which occur in this 
case. All the upper and lower case d-constants can be adjusted to 
satisfy the negative boundary conditions, but this will destroy the 
equilibrium of the central boundary solution since any change in the 
upper-case D's does effect the central equations. If one of the upper- 
case D's is constrained by the negative boundary, then the remaining 
unused lower-case c can be brought in to provide enough arbitrary 
constants to satisfy the central boundary conditions. The net effect of 
all this will generally be that a growing solution must be made to 
have a nonnegligible contribution at the central boundary. Conse- 
quently, it will be large at the positive boundary. This is illustrated 

in Fig. 6. 

When the situation discussed above and represented in Fig. 6 oc- 
curs, the residual values of y k outside the equalization region will be 
large and will generally grow larger as the number of taps is increased. 

As an actual example of a bad solution of the type discussed above, 
consider the pulse illustrated on page 563 of Ref. 5. The polynomial 



-Ic-i-fl 




Fig. 6 — Typical solution behavior when the conditions of equations (46) are 
not satisfied. 
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which represents this pulse with fS set as the peak value is 
B(z) = 15z 8 - 20z 7 + 40z 8 - 6O2 5 + Oz* + 115z 3 

+ 60z 2 + 20z + 10. (47) 
Notice that with this choice for fi , 



M = 5 m = 3. 



(48) 



The root locations for this polynomial are illustrated in Fig. 7. There 
are four inside roots (/ = 4) and four outside roots (O = 4). Con- 
sequently, equations (46) are not satisfied and no good truncated 
solution exists. This is verified by the fact that when Lucky at- 
tempted to equalize this particular pulse, the equalizer gave an output 
pulse with large y k ,k > N. This indicates that a solution such as that 
shown in Fig. 6 has been approached. 



i Im(z) 




Re(z) 



UNIT CIRCLE 



Fig. 7 — Root locations of pulse page 563 of Ref. 5. B(z) = 152 8 - 20z 7 + 40z e - 
60-z 8 + Oz* + 115z» + 60z* + 20z + 10 6 h ft* = 0.404 Z ± 89.8° 0», 6 2 * = 0.652 
Z ± 109.9°*!, *i* = 1-642 Z ± 27.7° 02, 02* = 1.894 Z ± 94.4°. 
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If p were shifted (delayed) one sample so that it would be the 
coefficient of the fourth-degree term of B(z), then equations (46) 
would be satisfied and a good truncated solution would exist. Notice 
however, that in this case /? would be zero which presents a severe 
problem in iteratively searching for this good solution. 

VI. ITERATIVE SEARCH FOR EQUALIZATION 

The discussion has indicated that, except for singular cases where 
one or more roots of B(z) are on the unit circle, equalization solutions 
exist (for the infinite tap case) in which tap weights decrease expo- 
nentially in both directions from the center tap. Furthermore, good 
truncated solutions exist which force all y k (— n ^ k ^ N, k 7± 0) 
to zero if ft is selected to satisfy equations (46), which can always 
be done. The question which now arises is whether a simple iterative 
search routine will lead to a desirable equalization. 

When presented with an impulse response to be equalized using a 
truncated equalizer and the Lucky criterion of forcing output samples 
to zero, a hierarchy of questions must be considered: 

(i) Is the pulse equalizable? That is, are all the roots off the unit 
circle? (If there are roots on or very close to the unit circle, a change 
of timing, that is, varying t in equation (3), will usually move the 
roots off the unit circle.) 

(ii) If the pulse is equalizable, does a good truncated solution exist? 
A shifting of (3 subscripts can always guarantee the existence of a good 
solution by satisfying equations (46) , but will sometimes create con- 
vergence problems. 

{Hi) If the pulse is equalizable and a good truncated solution exists, 
will a simple iterative search find this solution? 

The iterative method of searching for a solution which we consider 
first consists of measuring the value (y fc ) of the fcth output sample, 
then subtracting some part of this from the fcth tap weight.* That is 

«T» = «r - *rf r) (49) 

where the superscripts indicate the iteration number and A is a posi- 
tive number less than one. This iterative process is not identical to 
the method presented in Ref. 5 which increments the tap according 



♦This is not the only possible iterative search method, but it is one of the 
simplest. 
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to only the sign of y k . 



(r + l) . (r) 



a k 



= a y - Asm fori. 



(50) 



However, the two methods are very similar, and convergence of one 
will almost always guarantee convergence of the other. 

In matrix notation, the linear iterative search of equation (49) can 
be written as 



(r + l) (r) 



,(r) p (r) 

€ = na — n 



where 



«1 
LA-n J 



B = 



MO M-l M-2 

ft ft, /?-, /3_ 2 

ft ft ft /3_, /3_ 2 ••• 
•■• & /3, ft ^_! /3_ 2 
ft 



ft 








1 




0. 

o ••• o] 



0-, 



ft_ 



(51) 
(52) 



(53) 
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Equation (52) can be modified by using equation (51), 

e (r) - B[a (r ~ l) - Ac (, -°] - M 

= (7 - AJ5)e (r -" 

= (/ - A£)V 0> 

where I is the identity matrix. 
Clearly, the iterative routine will converge if 

lime (r) = 0. 



(54) 



(55) 



However, it does not quarantee convergence to a good solution as 
illustrated by the pulse of equation (47) and Fig. 7. There, con- 
vergence does occur as was illustrated in Ref. 5 and hence, equation 
(55) is satisfied but convergence is to a bad solution. In that case 
no good solution exists to which the routine can converge. 
The necessary and sufficient condition for convergence is that 



lim (I - AJ3) r = 0. 



(56) 



The matrix of equation (56) , that is, the matrix 7-AB will converge to 
zero if and only if all its eigenvalues are less than one in magnitude. 6 
This is equivalent to the condition, illustrated in Fig. 8, that the 
eigenvalues of B lie within a circle of radius 1/A centered at 1/A on 
the real axis. Assuming that A can be made as small as necessary, a 
necessary and sufficient condition for iterative convergence is that all 




REAL 
AXIS 



Fig. 8 — Necessary relation of A to location of eigenvalue in complex plane. 
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eigenvalues of B should have a positive real part. If this condition is 
not satisfied, {I— AB) r will diverge, and even if the initial error vector 
c (0) is very small, the final error vector will grow without bound. 

The italicized condition above, coupled with the root location con- 
dition of equation (46), guarantees convergence to a good solution. 
Certain other more restrictive sufficient (but not necessary) condi- 
tions can be derived. These conditions may be easier to check for a 
channel impulse response. 

6.1 Necessary and Sufficient Condition for Monotonic Convergence. 

The convergence of the iterative process will be called monotonic 
if for any starting error vector £ (0) , the following inequalities hold (a 
superscript t indicates the transpose of a matrix). 

He"" || > || i c " || > || e m || > ••• > \\e (r) || > \\e (r + l) \\ > ••• 

(57) 

|| f W || =e CD V „ = £ [6 (,> ]8 (5g) 

Since 

e (r + i, = ( 7 _ AB)e lr \ (59) 

it follows that 

|| e " +1 > || = e <-i>V r+1) = 6 (r, '(7 - AB)'(1 - AB)e (r) 

= e (r) V r) + A 2 [Be (r) ]'[Be (r) ] - Ae (r) '[B' + S]e (r) 

= || e (r) || + A 2 || Se (r> || - Ae (r, '[B' + B]e lr) . (60) 

Thus, equation (57) will be satisfied for all possible initial error vec- 
tors if an only if 

6*[£' + B]e > A || Be || = A(Be)'(Be). (61) 

Since the right side of the inequality is always positive, the inequality 
can be satisfied for all possible nonzero e only if [B* + B] is positive- 
definite. If this is true and 

then equation (61) will be satisfied for all * and the iterative process 
will be monotonically convergent. 
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If [B< + B] is not positive-definite, convergence can still occur. 
However, stating the conditions for convergence becomes more dif- 
ficult. For example, if \B f +B] is not positive-definite then for some 
vector, say € 0_1) , 



(r-l) l(nl 



[/?' + BW r ~" < 



and 



1 



II < 

\B\ j 

\B' -f 



(63) 

(64) 



+ 2A 2 || Be || 

£]e|| - 2Ae , [B' + B]e) (65) 



But if fore = £ ( '- n 

(Bt)'\B' + B\(Be) > ~ { [A 4 

+ A 1 
then 

|| e (r ~ l) || > || e (r + 1 || . (66) 

Very roughly speaking, equations such as (65) which can be developed 
indicate that [B'+B] should yield a predominantly positive quadratic 
form in order to have convergence. What we mean by this can best 
be illustrated graphically, as in Fig. 9. 

6.2 Further Sufficient Conditions jor Monotonia Convergence 

The positive-definiteness of [B'+B] is a necessary and sufficient 
condition for monotonic convergence. Being somewhat more restric- 
tive will yield other sufficient (but not necessary) conditions. We 
notice that [B l + B] is a Toeplitz matrix, that is, 



[B* + B] = [&„•], &„ = bu-.-i = 0„-„ + /3 ( ,-„ . 



(67) 




— NEGATIVE REGION 



VARIATION IN DIRECTION OF VECTOR e WITH 
ITS MAGNITUDE ||E|| HELD CONSTANT 



Fig. 9 — Illustrating :i predominately positive quadratic form. 
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Thus, it has some of the necessary attributes of a correlation matrix. 
A correlation matrix is positive-definite, so if we can show that 
[B'+B] is a correlation matrix, then this is sufficient for positive- 
definiteness. If [B' + B] is a correlation matrix then the polynomial 
B(z) + B l (z) = P (z), where 

B(z) = M z' v+m + • • • + /3_„, 
B\z) = (3. m z 2m + ■ • • + (3 Al z m - M (68) 

must have no odd order roots on the unit circle. Root locus methods 
may be applied here, that is, the loci of 



B(z) 
B\z) 



= -k 



(69) 



must not cross the unit circle. Notice that this is a necessary condi- 
tion that [B' + B] be a correlation matrix and it is sufficient but not 
necessary for [B' + B\ being positive-definite. 

A more restrictive condition leads to another sufficiency condition. 
If (Lucky condition, Ref. 5) 



0o > £ I 0* I- that is . D » < 1 



(70) 



then it is easy to see that B* (z) + B(z) has no roots on the unit 
circle and thus [B 1 + B] is a correlation matrix and positive-definite. 
Notice that equation (70) is not a necessary condition for [B l + B] 
being a correlation matrix. 
The various conditions discussed above are summarized in Fig. 10. 



.- SPACE OF ALL B MATRICES 




00 SATISFIES ROOT-LOCATION CRITERION, Eq.(46) 



-CONVERGENT REGION 
WHERE Re (X)>0 



MONOTONIC CONVERGENCE 
WHERE (Bt+B) IS 
POSITIVE - DEFINITE 

- (Bt + B) IS CORRELATION 
MATRIX 



"^—INITIAL ABSOLUTE DISTORTION D <1 

Fig. 10 — Illustrating the different convergence regions. 
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VII. MODIFIED ITERATION PROCESSES 



In the sample iteration process which has been discussed to this 
point, the change in tap weight «,• depends only on the output sam- 
ple value y;. Let us generalize so that the change depends on a linear 
combination of output samples. 



■ 7 



(r) (r> rj (r) 

e =7 — n — Ba — (J. 



e 
That is to say, the (r + list value of tap ; is 



(71) 



or, = a, -A2, v n .e k . (72) 

k 

Manipulation of equations (71) leads to a result similar to equation 
(54) 

€ (r) = [/ - AB7]V 0) . (73) 

Consequently, all the sufficient and necessary conditions which have 
been developed on the previous pages can now be applied to the 
matrix BV. Now, however, we have considerably more latitude since 
we are free to specify V. 

As an example, let us suppose that V is chosen to equal B*. Then 

BV = BB l = positive-definite (74) 

and monotonic convergence is guaranteed. This is a particularly ap- 
pealing way of selecting V since the sample weighting can be deter- 
mined directly from the initial channel impulse response. 

v ik = » f _ t = &_,- . (75) 

It is very interesting that the weighting suggested above is very 
nearly equivalent to inserting a tapped-delay-line matched filter 
ahead of the equalizer. A matched filter, whose tap weights are equal 
to the ft values in reverse order, will yield an output whose samples 
will form a B matrix which is a correlation matrix. Thus, the itera- 
tive search will be monotonically convergent in this case also. The 
weighting suggested above yields the same matrix except for "edge 
effects." 

This can be illustrated by the following example. Suppose there 
are just three /? values: /?_i =1, /? = 2, fii = —2. Then, with V set 
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equal to B f , and assuming a six-tap equalizer so that B is a 6 X 6 
matrix, we obtain 







-2 

-2 

8J 

On the other hand, the modified B matrix following a matched 
filter would be 



BV = BB l = 



5 


— 2 


-2 








-2 


9 


-2 


-2 





— 2 


-2 


9 


-2 


-2 





— 2 


-2 


9 


-2 








_2 


_2 


9 











2 


-2 



(76) 



tf, n = 



9 


-2 


—2 











— 2 


9 


-2 


-2 








-2 


-2 


9 


-2 


— 2 








— 2 


— 2 


9 


— 2 


-2 








_2 


— 2 


9 


— 2 











_2 


_2 


9 



(77) 



Setting V — B t will guarantee monotonic convergence. However, 
since in most cases B is such that [B l + B] is close to being positive- 
definite, it is probable that a less extensive V would be sufficient to 
guarantee monotonic convergence. As an example, suppose that V 
is chosen to be a small deviation on the standard iteration of equa- 
tion (58) 



7 = 7+ SB 1 , 



8 > 0. 



Now, 



BV = B + 8BB' 



(78) 



(79) 



will have a quadratic form which is greater for every vector c than 
the quadratic form for the matrix B alone. Thus, if the negative re- 
gion such as is illustrated in Fig. 8 is small, then BV can become posi- 
tive-definite for relatively small S. 

Perhaps a more reasonable way of selecting a V which approximates 
B t is to modify equation (75) in the following manner. 
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.„-*-- (*- if , ^'- L (80) 

[0 if I p k - t | < L 

Thus, only the more significant values of /? are used in weighting the 
errors. In general such a V will force B towards positive-definiteness. 
However, there seems to be no general rule for selecting the critical 
value of L, that is the maximum L which will just permit positive- 
definiteness. 

VIII. CONCLUDING REMARKS 

Virtually all input pulses are equalizable in the sense that there 
exist tap-weight adjustments which will force the output samples in 
the adjustment interval to zero while the output samples out of the 
adjustment interval remain small. Furthermore, the residual samples 
outside of the adjustment interval will become smaller as the number 
of taps (length of transversal filter) is increased. 

There are just two necessary conditions in order for the preceding 
statement to hold. The first is that the polynomial representing the 
input pulse have no roots on the unit circle. Although the singular 
case where roots are exactly on the unit circle is highly improbable, 
roots very near the unit circle lead to relatively larger residual er- 
rors and greater potential for instability. 

The second necessary condition is that the selection of the central 
sample value must be such that equation (46) is satisfied. 

Although a pulse may be equalizable, the simple first-order itera- 
tive search for the proper tap weights given by equations (51) and 
(52) may not be convergent, If it is convergent, and assuming the 
two conditions above are satisfied, it will converge to the proper tap- 
weight settings. If it is not convergent, it will be divergent with in- 
creasing errors in the adjustment interval. The convergence or diver- 
gence is independent of the initial tap settings. Thus, even though 
the tap weights might be set to optimum initially, if the system is in 
the iterative search mode and is divergent, it will eventually diverge. 

The necessary and sufficient condition for convergence is given in 
Section VI along with a hierarchy of more stringent sufficient condi- 
tions. In general, convergence will be dependent upon the absolute 
timing of the sampling. Consequently, a particular pulse which is 
equalizable for two different timings might be convergent for one 
timing and divergent for the other. 

If the first-order iterative procedure is divergent, a more complex 
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weighting of the output errors in adjusting the tap weights can im- 
prove the situation. At least one weighting given by equation (75) 
will guarantee convergence. 
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